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cN : Abstract 

^ ! Explaining the formation and evolution of galaxies is one of the most 

^ ■ challenging problems in observational cosmology. Many observations sug- 
gest that galaxies we see today could have evolved from the merging of 

0^ ■ smaller subsystems. Evolution of galaxies tells us how the mass or number 
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density of the lens varies with cosmic time. Merging between the galaxies 
and the infall of surrounding mass into galaxies are two possible processes 
that can change the comoving number density of galaxies and/or their 
■ mass. We consider five different evolutionary models of galaxies .These 

models are: Non evolutionary model, Guiderdoni and Volmerange model, 
fast merging, slow merging and mass accretion model. We study the grav- 
itational lens image separation distribution function for these models of 
evolving galaxies. A comparison with data for lensed quasars taken from 
the HST Snapshot Survey rules out the fast merging model completely 
as this model produces a large number of small-separation lenses. It is 
possible that the mass accretion model and the non evolutionary model of 
galaxies may be able to explain the small angle separations. 
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1 Introduction 



After the discovery of the first multiply imaged quasars, gravitational arcs 
and arc-lets, gravitational lensing has rapidly become one of the most 
promising tools for cosmology. It is not a new idea that the statistics 
of gravitational lensing can be used as a tool for the determination of 
cosmological parametersJ'S In their pioneering work. Turner, Ostriker 
& Gott (1984, hereafter TOG)i developed a formalism to calculate the 
lensing probability and image separation distributions. TOG modelled the 
lens population as point masses, singular isothermal sphere (SIS) galaxies 
and cluster of galaxies. Since gravitational lens frequencies are sensitive 
to the cosmological constant (A), many author^^ used the lens statistics, 
developed by TOG, to constrain the cosmological parameters. Hinshaw 
& KraussS (1987) and Krauss k White0(1992) studied these effects with 
finite core radii of galaxies. In all these papers one simplifying assumption 
is made that the comoving number density of galaxies (lenses) is constant. 
However, it is an oversimplification to assume that galaxies are formed at 
a single epoch. 

Maolilf 1991) first examined the effect of galaxy evolution on the statis- 
tical properties of gravitational lenses using a simple redshift cut-off model. 
This model reduces the lensing probability and also explains the large 
mean separation of images in the observed gravitational lenses. Sasaki & 
Takaharali(1993) used a more realistic model than that of Mao and stud- 
ied the effect of a more gradual evolution in the number of galaxies on lens 
statistics. Mao & KochanekS (1994) put a hmit on galaxy evolution by 
studying its effects on gravitational lens statistics and image separations. 
They concluded that most of the galaxies must have collapsed and formed 
by z ~ 0.8 if the universe is described by the Einstein-de Sitter model. 
If elliptical galaxies are assembled from merging of spirals then most of 
ellipticals must be formed by a redshift of 0.4. 

Rix et al. 1(1994) considered an evolutionary model which is physically 
more plausible. They have studied the gravitational lens statistics with 
an evolving lens population through the fast merging process only. They 
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showed that some specific merger models can be rejected and A- dominated 
cosmologies are ruled out as they predict a large number of sub-arc second 
lenses and with merging, this problem becomes more acute. They find that 
the lensing probability of getting multiple images is insensitive to merging 
and the merging scenario skews dN/d{A6) towards smaller separation. 
Park & Gottlil(1997) tried to explain the correlation between gravitational 
lens image separation and the source redshift in the presence of the galaxy 
evolution. Jain et al.0 (1998a) also studied the effect of galaxy evolution 
on the statistical properties of lenses both with decaying and constant A. 

The aim of this paper is to use the image separation distribution func- 
tion of lensed quasars, {dN/dA9), as a tool to put constraints on the 
various evolutionary models of galaxies. Hamana et al. ^(1997) studied 
the distribution of image separation angle of lensed quasars. However 
they found a complete mismatch between the theoretical estimates and 
the observed lensed events. We find that if we use {dN/dA6) as a tool to 
study galaxy evolution, it not only reduces the gap between observations 
and theoretical predictions but also tells us which evolutionary model of 
galaxies might explain the observations. The paper is organized as fol- 
lows. In §2 we explain the different evolutionary models of galaxies. In 
§3 we present the statistical formulas we use for comparing the model 
calculations with the observations. We briefly summarize our results in 
§4. 

2 Evolution Of Galaxies 

The theory of the formation and evolution of galaxies is one of the un- 
solved problems of astrophysics. Some authors believe that galaxies evolve 
through a complex series of interactions before settling in the present day 
form.lll'll^ Others believe that galaxies were created in a well defined event 
at very early time.iil'il It remains unclear which process dominates the 
formation of elliptical galaxies. Among the many theories of galaxy for- 
mation, the idea that galaxies may form by the accumulation of smaller 
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star forming subsystems has recently received much attention. Many ob- 
servations also support this 'bottom-up' scheme. 

First, deep Hubble Space Telescope (HST) images^l indicate that early 
type galaxies were assembled largely at 2; > 1 and have been evolving 
passively since z < 1. Moreover, HST and ground based telescopes show 
that the galaxy-merger rate was higher in the past and it roughly increases 
with redshift.S^ This suggests that the galaxies we see today could 
have been assembled from the merging of smaller systems sometime before 
z ~ 1. 

Recent observations^^ also show that elliptical galaxies are rarer at high 
redshifts than those predicted by models in which elliptical galaxies com- 
pleted their star formation by 2; ~ 5. Therefore elliptical galaxies must 
have had significant star formation at z < 5 through merging and asso- 
ciated starbursts. The formation of elliptical galaxies in this way is also 
consisitent with the predictions of hierarchical clustering models of galaxy 
formation. 

The second piece of evidence comes from the excess of faint blue galax- 
ies which has been found in many deep imaging studies.^3^11 Comparison 
with the model which assumes that no luminosity evolution takes place 
in the galaxy population, shows that in the B band the actual observed 
galaxy count exceeds the model predictions by a factor of 5. Merging of 
galaxies can solve the surprisingly steep increase in the number density 
of galaxies. i^'il But at present it is not clear whether such models ade- 
quately describe the merging of galaxies in any realistic models of structure 
formation in the universe. 

We consider five different evolutionary models of galaxies which try to 
explain some of the observational facts listed above . These models are: 
Non evolutionary model, Volmerange and Guiderdoni model, fast merging, 
slow merging and mass accretion model. We study the gravitational lens 
image separation distribution function in the presence of evolving models 
of galaxies. 
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Non- Evolutionary Model 



This is the conventional lens model in which the luminosity function of 
lens galaxies is assumed to be of the Schechter form^^ 

$(L, ^ = 0) = (J)*{L/L,)''exp{-L/L,)dL/L, (1) 

where the 0*, a and L* are the normalization factor, the index of faint 
- end slope, and the characteristic luminosity respectively. These values 
are fixed in order to fit the current luminosities and densities of galaxies. 
This model assumes that the comoving number density of galaxies n{6t) 
is constant and the mass of galaxies does not change with cosmic time. 

n{St) = no = constant (2) 

where 6t is the look-back time. The velocity dispersion of Singular Isother- 
mal Sphere (SIS) lenses at 6t is 

v{6t) = vq = constant (3) 
The subscript refers to present-day values. 

Volmerange and Guiderdoni Model 

In 1990, Volmerange and Guiderdoni,^^ proposed a unifying model to 
explain faint galaxy counts as well as observational properties of distant 
radio galaxies. This model of galaxy evolution is based on number evolu- 
tion in addition to pure luminosity evolution. According to this model the 
present day galaxies result from the merging of a large number of build- 
ing blocks and the comoving number of these building blocks evolves as 

It is argued that the present luminosity function is the well known 
Schecter Luminosity Function^ given in eq.(l) above. Then at high z, the 
comoving number density follows New Luminosity Function 

$(L, z)dL = (1 + z)2'^$(L(l + zy, 0)dL (4) 
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It is seen that the value 77 = 1.5 gives a fair fit to the data on high redshift 
galaxies. The functional form has the following properties: 

(i) Self-similarity as suggested by the Press-Schecter (1974)111 formalism 
subject to the constraint that the total mass of associated material is 
conserved. 

(ii) The comoving number density evolves as (j)*{z)=(j)Q (1 + z)^ and the 
characterstic luminosity of the self similar galaxy luminosity function 
varies as L^(z) = L^(0)(1 + z)^^ 



Fast Merging Model 

The first merger model is that of Broadhurst, Ellis & Glazebrook (1992),!: 
which was originally motivated by the faint galaxy population counts. This 
model assumes the number density of the lenses to be a function of the 
look back time 6t as: 



n{5t) = f{5t)no (5) 
The velocity dispersion of SIS lenses at 6t is 

v{6t) = [f{St)]-\ (6) 

This form implies that if we had n galaxies at time 6t each with velocity 
dispersion they would by today have merged into one galaxy with a 
velocity dispersion [f{6t)]v. The strength and the time dependence of 
merging is described by the function f{St): 

f{St) = exp{QHo6t) (7) 

where Hq is the Hubble constant at the present epoch and Q represents 
the merging rate. We take Q = 4.il The look back time 6t is related to 
the redshift z through 
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where F{y) = Qo(l + + (1 - - Ao)(l + y)^ + Ao , 

Qq = 87tGpq/3Hq , Pq is the density of matter, Aq = SirGpyo/SHQ and pyo 

is density of vacuum energy. 

Slow Merging Model 

In this less extreme merger model the total mass of galaxies within a given 
comoving volume is conserved. The comoving number density goes like 
while the mass of an individual galaxy increases like t^/^, where t 
is the cosmic time since the big bang.li^'^^ We further assume the mass- 
velocity relation M (x v'^. The values of 7 and are given in Table 2 for 
elliptical galaxies. Then 



n{6t) = no 



1-^ 



v{St) = Vq 
where to is present age of the universe. 



1-^ 



-2/3 



37 



(9) 



(10) 



Mass Accretion Model 



Mass accretion is the key factor for evolution of galaxies. A galaxy can ac- 
crete mass through two processes: either it accretes gas regularly through 
internal dynamics or the accretion occurs in more violent events, galaxy 
interactions and mergers. This is in line with the idea of hierachical for- 
mation. In this model the comoving density of the galaxies is constant but 
the mass increases as 

^2/3 

as in the cosmological infall model.0 The total 
mass in galaxies thus increases with time. The comoving number density 
and the dispersion velocity vary as 



n{6t) = no{constant) 



11: 



v{St) = Vq 



6t- 



37 



(12) 
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3 Basic Equations For Gravitational Lensing Statis- 
tics 



The differential probability dr of a beam encountering a lens in traversing 
the path of dzL is given by 

cdt 

dr = nL{z)a- — dzi^ (13) 

dZL 

where til^z) is the comoving number density [TOG 1984]. 

The Singular Isothermal Sphere (SIS) provides us with a reasonable 
approximation to account for the lensing properties of a real galaxy. The 
lens model is characterized by the one dimensional velocity dispersion v. 
The deflection angle for all impact parameters is given hy a = Attv'^/c^. 
The lens produces two images if the angular position of the source is less 
than the critical angle Per: which is the deflection of a beam passing at 
any radius through an SIS: 

(3,r = aDLs/Dos. (14) 

we use the notation Dql = d{0, zl), D^s = d{zL, zs), Dqs = d{0,zs), 
where d{zi^ Z2) is the angular diameter distance between the redshift zi 
and Z2^ Then the critical impact parameter is defined by = DoiPcr 
and the cross- section is given by 

-^'^'-''^'&(^)^' (^^) 

The Evolutionary Model 



The differential probability dr of a lensing event in an evolutionary 
model can be written asi^'il 
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dr = 



DnrD 



OLJ^LS 



1 CC?f 



dZL 



(16) 



Ro dzL 



dzi 



where F = ^^<i>*vjr (a + ^ + 1^. The functional form of f{St) in 
various models is described in Table 1. 



Table 1 

The Functional Form of f(5t) 



Evolutionary Model 


Form of (fSt) 


Fast Merging 
VG Model 

Slow Merging 
Mass Accretion'' 


exp{QHo5t) 
(1 + ^)1-5 

(1-1)"" 



superscript 'a': In this case the exponent of f{5t) in eq. (16) becomes (—1 — ^) as the 
total mass in galaxies increases with time. 

Using the relations = g^f^^^, = 87r(^)^^^ and luminosity 

- velocity relation = {^y, we get = (^^0) 2 . is image angular 
separation. The differential optical depth of lensing in traversing dzL with 
angular separation between and d(j) is given by 



dzidcf) 



d(t)dzL = F{1 + z)^exp 



Dos 
Dls 



X 



DosRo 
Dos 
Dls 



Ro dzL 
( 



7/2 



—dzj 



(17) 
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In eq.(17) the exponent of f{St) becomes a for the mass accretion model 
while in other evolutionary models the exponent of f{St) remains a + 2. 
After integrating eq.(17) over the lens redshift zl, we obtain the angular 
separation distribution 



dN 



=<5>E 



d^r 



-dzi 



(18) 



dAO ^Jo dzLdAe 

where the summation is over sources in the HST quasar sample and 
< B > is the averaged bias which is equal to 9.76 in the HST quasar 
sample.El 



The Non Evolutionary Model 

In the non- merging model the optical depth is given by! 



dr = 



vtr (a + i + l) (1 + zl] 



DnrD 



X 



OLJ^LS 



4 

7 

,^ 1 cdt 



dzL 



(19) 



RoDqs / dzL 

It is clear from eq. (16) that if 7 = 4 then optical depth at each redshift 
in evolutionary model of galaxies (except mass accretion model) becomes 
equal to the optical depth in the non evolving model of galaxies. The 
differential optical depth of lensing in traversing dzi for angular separation 
between and (j) + d(j) for the non evolutionary model is given by 



d\ 
dzidcj) 



d(t)dzL = F{1 + z)^exp 



D 



OS |^2 



Dls 



x( 



X 



DolDls ^2 

DqsRo 

Dos 



1 cdt 
Ro dzL 



Dls 



T(a+ 1 + 



—dzi 



(20) 



The dN/d{A6) which is given by eq.(18), strongly depends upon the 
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four parameters a, 7, 0* and dispersion velocity v^. We consider three 
sets of these parameters given in Table 2. 



Table 2 



Schechter and Lens Parameters for E/SO Galaxies 



Survey 


a 7 v*{Km/s) (j)*{Mpc~^) F* 


K96 

LPEM 

NS97 


-1.0 4.0 225 6.1 X 10-3 0.026 
+0.2 4.0 205 3.2 x 10"^ 0.010 
-1.1 3.3 175 1.1 X 10-2 0.015 



References: K96 - C. S. Kochanek (1996).! LPEM - J. Loveday et al. (1992). i NS97 - 
T.T. Nakamura & Y. Suto (1997)@ 



4 Result and Discussion 

An obvious noticeable fact from eq. (16) is that when the index 7 (the 
Faber- Jackson index) is different from the value of four that we get any de- 
pendance on the evolution of galaxies. The combination (1 — 4/7) vanishes 
when 7 has the value 4 and totally suppresses the effect of evolution. 

Figs. 1 to 3 give the calculated and observed number of lensed objects 
as a function of the image separation for the currently accepted values of 
cosmological parametersSil Q = 0.3 and A = 0.7. In all these figures the 
histograms indicate the image separation distribution of the four lensed 
quasars observed in the HST snapshot survey. It is clear from these figures 
that there is a very significant dependance of the results on the Schecter 
and lens parameters. The uncertainties associated with these parameter^ 
will be diminished greatly by the data that emerges from the surveys 
planned for the next few years. The K96 set of parameters shift the peak of 
the image separation distribution function ( dN / d{A6)) for all the evolving 
models towards smaller angles of separations. The larger angle separations 
are not explained by these parameters. With LPEM parameters, the mass 
accretion model and, to a lesser extent, the non-evolution model seem to 
have reasonable success in explaining the observations. Figs 4 to 6 present 
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the results for the case of zero cosmological constant. Here the best fit 
seems to be the K96 parameters. 

It is clear that to reach at a firm conclusion about galaxy evolution 
more reliable and valid set of Schechter parameters are badly needed. 
This should be forthcoming in the next few years with the large number 
of surveys that are in progress and that are being planned for the future. 

Satisfactory and unambiguous identification of lensed objects should 
also improve in the coming years. The X-ray observations by CXO (Chan- 
dra X-ray Observatory) should be specially helpful.@ An increase in the 
overall number of lensed objects at all wave lengths is also expected in 
the next few years. It is hoped that knowledge of cosmological parameters 
and galaxy distribution (Schechter) parameters will mutually refine each 
other as more results become available. It will then be possible to say with 
some certainty which of the various galaxy evolution models are favoured 
by observations and which are ruled out. 
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Q = 0.3, A = 0.7 
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Figure 1: 

expected distribution of lens image separations with K96 parameters with Q 
A = 0.7 
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a = 0.3, A = 0.7 




Figure 2: 

expected distribution of lens image separations with LPEM parameters with Q 
A = 0.7 
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Figure 3: 

expected distribution of lens image separations with NS97 parameters with Q 
A = 0.7 
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£2 = 1 , A = 

T — I — I — r 
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Figure 4: 

expected distribution of lens image separations with K96 parameters with Q 
A = 0.0 
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£2 = 1 , A = 



■ ■ Mass Accretion 

VG Model 

Fast Merging 

Slow Merging 

No Evo. 




Figure 5: 

expected distribution of lens image separations with LPEM parameters with Q 
A = 0.0 
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£2 = 1 , A = 
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Figure 6: 

The expected distribution of lens image separations with NS97 parameters with Q 
1.0, A = 0.0 
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